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f^) , The energy spectra of a variety of collective modes on quasiperi- 

' odic structures exhibit a complex fractal profile. Among the modes 

that have attracted particular attention in this context, are the spin 
wave spectra of quasiperiodic magnetic multilayers that obey a sub- 
stitutional sequence of the Fibonacci type. They are described within 
' the framework of the Heisenberg theory. In order to have a deep in- 

sight on the relevant thermodynamical implications of the above men- 
tioned energy spectra's fractal profile, we have performed analytical 
and numerical calculations of the spin wave specific heat associated 
with successive hierarchical sequences of the Fibonacci quasiperiodic 
, structures. The spectra show interesting oscillatory behavior in the 

low-temperature region, which can be traced back to the spin wave's 
self-similar energy spectrum. 
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1 Introduction 



Since the discovery of the icosahedral phase in Al-Mn alloys by means of 
X-ray spectroscopy the quasicrystaline systems have been extensively 
studied (for a revision see [§]). In particular, the physical properties of the 
so-called quasiperiodic structures have attracted a lot of attention from both 
theoretical and experimental point of view. Although the terni quasicrystal 
is more appropriate for naturai compounds or artificial alloys, in one dimen- 
sion there is no difference between this case and the quasiperiodic structure 
formed by the incommensurate arrangement of periodic unit cells. Due to 
this motivation, Merlin and collaboratore, using the molecular beam epitaxy 
(MBE) technique, grew in 1985 the first quasiperiodic superlattice following 
the Fibonacci sequence H. After that pioneer work, other quasiperiodic 
structures were experimentally realized ||]. 

^From a theoretical point of view, a number of physical properties has 
been studied in quasiperiodic structures. Among them we can cite the energy 
spectra of polaritons || , phonons , electrons Q and spin waves || |9| , as 
well as the magnetoresistance and magnetization curves of quasiperiodic thin 
films [10]. A quite interesting feature, common to ali of these systems, is a 
self-similar pattern of their spectra. In fact, the energy spectra of the above 
referred particles and systems are highly fragmented and tend to Cantor sets 
in the thermodynamic limit. The origin of this fractality can be attributed 
to the long range order induced by the non-usual hierarchical structure of 
the quasiperiodic sequences used in the construction of the system. 

A very interesting question aroused in the last few years: what are the 
consequences of a fractal energy spectra on the behavior of the thermody- 
namic properties of quasiperiodic systems? In a recent work, Tsallis and 
collaboratore [[il]] studied the specific heat properties of a fractal energy 
spectra generated by the geometrical triadic Cantor set. They have found, 
as the main result of the paper, that the specific heat presents oscillations 
around the fractal dimension of the spectra. In addition, a non-uniform con- 
vergence between the so-called banded and discrete models was observed. 
Their results were extended by Vallejos et al [O, [l3j for the two-scale Can- 
tor set case. For this more general situation, the specific heat also exhibits 
log-periodic oscillations around the fractal dimension. Later Carpena et al 
[14|, using the properties of a multifractal spectra showed under what condi- 
tions the oscillatory regime disappears. Finally, Curado and Rego-Monteiro 
[15| examined the thermodynamic properties of a solid exhibiting the energy 
spectrum given by a logistic map. 
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The aim of this work is to push a little bit more the understanding of 
the thermodynamic properties of quasicrystals, analyzing the specific heat 
of a real fractal spectra of spiri waves in quasiperiodic Fibonacci magnetic 
superlattices. We consider that the building blocks used here to set up the 
quasiperiodic structures are ferromagnetic materials, whose dynamics are 
described by the Heisenberg Hamiltonian. Our main concern is to emphasize 
the differences between the specific heat profile obtained by using a real spin 
wave's multifractal energy spectra for a given in-plane wavevector, (with, by 
instance, their proper scaling laws), as discussed in Ref. 8, and an idealized 
Cantor set, considered in previous works. Throughout the paper we use the 
classical Maxwell-Boltzmann statistics. 

This paper is organized as follows: In section II we discuss the theoret- 
ical method used to obtain the spin waves dispersion relation (their energy 
spectra), and the rules used to build up the Fibonacci super lattice. Sec- 
tion III is devoted to the specific heat theory of a general scaled energy 
band spectra. The application of this theory to the Fibonacci quasiperiodic 
structure is then made, and the main features of the spectra obtained are 
discussed. Finally, section IV presents the conclusions of this work. 

2 Physical model 

In this section we briefly describe the physical model for the quasiperiodic 
magnetic superlattices. A detailed description can be found in Ref. ||. We 
consider superlattices in which uà layers of material A (building block A) al- 
ternate with iiB layers of material B (building block B). Both materials are 
taken to be simple cubie spin-S 1 Heisenberg ferromagnets having exchange 
constants J a an d Jb, respectively, and lattice Constant a. At the interface 
A/B, the exchange Constant is suppose to be equal to /. The Heisenberg 
Hamiltonian for each component is: 



Here the sum in the first term is over sites i and nearest neighbors (n.n.) 
j, and Hq is a static external magnetic field pointing in the z-direction . 

The spin wave dispersion relation in a superlattice can be found by solv- 
ing the equations of motion for the operator = Sf + iSf , i.e: 
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h-Sf = g^ B H Q SÌ+ <S Z >Y J JA{Sf - S+), (2) 

n.n. 

where < S z > is the random phase approximation for the z- component of 
the spin operator. The solution of (2), for material A, is: 

= [Ai exp(i kA • ~ta) + Bi exp(-i kA • r\4)](exp -io;*) (3) 

with similar expression for material B. These solutions are linked together 
using the equation of motion (2) at the boundaries of the n-th unit celi, 
which can be written in a matricial form. Using the translational invari- 
ance property of the excitations, through Bloch's theorem, the spin wave 
dispersion relation follows as: 

cos(QD) = (l/2)Tr [T] . (4) 

Here Q is the Bloch'wavevector of the collective mode, D is the size of 
the superlattice's unit celi, and T is a transfer matrix which relates the 
coemcients of the (l + l)th celi to the coemcients of the preceding one. 

A Fibonacci superlattice can be grown experimentally by juxtaposing 
the two building blocks A and B in such way that the ni/i-generation of the 
superlattice S n is given iteratively by the rule S n = S n -iS n -2, for n > 2, 
with Sq = B and Si = A. It is also invariant under the transformations 
A — > AB and B — > A. The Fibonacci generations are, 

So = [B], S l = [A], S 2 = [AB], S 3 = [AB A], etc. (5) 

The number of building blocks increases according to the Fibonacci number, 
F n = -F n _i + F n _2 (with i<o = -Fi = 1), and the ratio between the number of 
the building blocks A and the number of building blocks B, when n » 1, 
in the sequence tends to r = (1/2) (1 + \/5), an irrational number known as 
the golden mean. 

It can be shown that the transfer matrix for the nth generation of a 
Fibonacci superlattice can be obtained by a simple recurrence relation given 

by§, 

T an =T an _ a -T 8n _ 1 , n>2. (6) 

Therefore, from the knowledge of the transfer matrices T$ and Ts x , we 
can determine the transfer matrix for any generation and consequently the 
dispersion relation. Note that the matrix T$ 2 recovers the periodic case. 
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3 Specifìc heat Spectra 



The spin wave fractal spectra for the Fibonacci superlattices is depicted in 
Fig. 1, for a fixed value of the in-plane dimensionless wavevector k x a. From 
there we can see the forbidden and allowed energies of the spin wave spectra 
against the Fibonacci's generation number n, up to their 8th generation, 
which corresponds an unit celi with 21 ^4's and 13 i?'s building blocks. The 
number of allowed bands is equal to three times the Fibonacci number F n , 
of the correspondent generation. Notice that, as expected, for large n the 
allowed band regions get narrower and narrower and they have a typical 
Cantor set structure. 

We address now to the specifìc heat of the spectra indicated in Fig. 1. 
The description below is general and can be applied to any banded spectrum. 
In Fig. 1 each spectrum, for a fixed generation number n, has m allowed 
continuous bands. We consider, without loss of generality, the level density 
within each band Constant. The partition function for the nth generation is 
then given by: 



Here (5 = l/T (by choosing the Boltzmann's Constant ks = 1), and we take 
the density of states p(e) = 1. After a straightforward calculation we can 
write Z n as, 



Here the subscript n is the generation number, m is the number of allowed 
bands and Aj = e^+i — ej is the difference between the top and bottom 
energy levels of each band. 

The specifìc heat is then given by, 




(7) 



-, 2m-l 



i=l,3,... 



■13 Ai 



(8) 



C„(T) 



d 2 <5 m Zi 



n 



(9) 



dT l dT 



which can be written as, 



C n {T) 




(10) 
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with, 

2m-l 

fn= E k 2 e-^- e f +1 e-^]. (11) 

i=l,3,... 

and, 

2m-l 

<7n = E [e,e-^-e m e-^+ 1 ]. (12) 

i=l,3,... 

Therefore, once we know the energy spectra of the spin wave which 
propagates in a given sequence's generation of a quasiperiodic structure, we 
can determine the associated specific heat's spectra by using (10). 

Fig. 2 shows the spin wave specific heat spectra of the Fibonacci su- 
per lattices, for the in-plane wavevector k x a = 2.0, as a function of the 
temperature. For the high temperature limit (T — ► oo), the specific heat 
for ali generation numbers converges and decays as T -2 , for arbitrary n, in 
agreement with the triadic case. This is a consequence of the existence of 
a maximum energy value in the spectrum (once the spectrum is bounded). 
As the temperature decreases, the specific heat increases up to a maximum 
value. The corresponding temperature for this maximum value depends on 
the Fibonacci generation number n, although one can see a clear tendency 
for a common temperature value as n increases. After the maximum value, 
the specific heat falls into the low temperature region. In this region it starts 
to present non-harmonic small oscillation behavior, as shown in the inset of 
Fig. 2. These oscillations, are not around a specific fractal dimension of 
the quasiperiodic structure, as in the idealized triadic Cantor set! Besides, 
it cannot be considered also an approximation of the idealized oscillations 
found in the triadic Cantor set. Their profiles define clear ly two classes of 
oscillations, one for the even and the other for the odd generation numbers 
of the sequence, the amplitude of the odd oscillations being bigger than the 
amplitude of the even one. These behaviors are better illustrated in Fig. 
3, where are depicted log-log plots of the specific heat against the temper- 
ature for several generation numbers. Of course the number of oscillations 
observed in the specific heat spectra is related to the hierarchical generation 
number n (more oscillations appear as n increases). Another behavior of 
these oscillations can be seen in Fig. 4, where we have plotted the specific 
heat against logT. We can note a well defined period in the oscillations, 
which means that the specific heat is a log-periodic function of the temper- 
ature. The curves resemble the triadic case, with a mean value d, around it 
C(T) oscillates log-periodically, although, as in the log-log plot, this value is 
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not related with the fractal dimension of the Fibonacci quasiperiodic struc- 
ture. Besides, the mean value is different for the even and odd Fibonacci's 
generation numbers. 



4 Conclusions 

In this paper we have studied the spin wave sperine heat, for fixed values 
of the in-plane dimensionless wave vector k x a, for the Fibonacci magnetic 
superlattices. We have shown, as a common aspect of the model employed 
here, that the sperine heat tends to zero in the high temperature limit (T — * 
oo), as T~ 2 , no matter its generation number. This asymptotic behavior is 
mainly due to the fact that we have considered our system bounded. It 
also presents a maximum whose corresponding temperature tends to a fixed 
value as the generation number increases. Finally, in the low temperature 
region, a small oscillations arises as an indication that the sperine heat is a 
log-periodic function of the temperature. These oscillations can be defined 
as the signature of the quasiperiodic system, and has no counterpart in the 
idealized case. 

It would be of interest to have experimental data to test our predicted 
theoretical results presented here. However, most of the experimental stud- 
ies of magnetic multilayers to date, have been done for transition metals 



[16, 17]. Further experimental studies carried out for multilayers structures 
of magnetic insulator or magnetic semiconductor materials, that would be 
better described by the Heisenberg model employed in this paper, are wel- 
come. Suitable experimental technique to probe the spin wave spectra in 
the quasi-periodic structure discussed here is the inelastic light scattering 
spectroscopy of Raman and Brillouin type. However, the spectra can be 
obtained only for a given in-plane wavevector k x , which defines the incident 
angle 6 of the light through the relation sin -1 9 = k x \/4ir, A being the the 



laser wavelength [18, [19] . This fact adds an extra complication to the exper- 
imental specific heat measurements, but we hope that the experimentalists 
can be encouraged to overcome it. Techniques involving magnetic resonance 
(for example, ferromagnetic resonance, standing spin-wave resonant etc.) 
can also be used, and indeed they were previously been successfully applied 
to surface and bulk spin waves in various magnetic microstructures (for a 
good account of these techniques see f20|). 
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Figure Captions 

1. The spin wave spectra for the Fibonacci structure. Here the in-plane 
wavevector is taken to be equal to 0.2. 

2. Sperine heat versus temperature for the Fibonacci structure. We have 
plotted the sperine heat profiles up to the 9th generation number. The 
inset shows the low temperature behavior of the specific heat. 

3. Logdog plot of the specific heat versus temperature for the generation 
numbers of the Fibonacci quasiperiodic sequence. Observe a different 
behavior for the even (n=4, 6, 8, and 10) and odd (n=5, 7, and 9) 
generation numbers. 

4. C(T) vs log(T) plot to show the log-periodicity effeets in the Fibonacci 
structure. 
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